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ABSTRACT 

We present a time series of synoptic images of the hnearly-polarized v = 
1, J = 1 — SiO maser emission toward the Mira variable, TX Cam. These data 
comprise 43 individual epochs at an approximate biweekly sampling over an op- 
tical pulsation phase range of (p = 0.68 to (f) = 1.82. The images have an angular 
resolution of ~ 500/ias and were obtained using the Very Long Baseline Array 
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(VLBA), operating in the 43 GHz band in spectral-line, polarization mode. We 
have previously published the total intensity time series for this pulsation phase 
range; this paper serves to present the linearly-polarized image sequence and an 
associated animation representing the evolution of the linear polarization mor- 
phology over time. We find a predominantly tangential polarization morphology, 
a high degree of persistence in linear polarization properties over individual com- 
ponent lifetimes, and stronger linear polarization in the inner projected shell 
than at larger projected shell radii. We present an initial polarization proper 
motion analysis examining the possible dynamical influence of magnetic fields in 
component motions in the extended atmospheres of late-type, evolved stars. 

Subject headings: masers - polarization - stars: magnetic fields - stars: individual 
(TX Cam) 



Introduction 



The extended atmospheres of late-type, evolved stars on the asymptotic giant branch 
(AGB) are host to complex underlying astrophysical processes of vital importance. Mass- 
loss from these stars originates in this region and is an important enrichment mechanism 
for the interstellar medium. The large- amplitude, long-period variables (LALPV) span a 
range of variability sub-cl asses, including the Mira variables, with a median pulsation period 
of several hundred days (IHabing I Il996l ) . Their near-circu mstellar environments (NCSE) 
are permeated by pulsation shocks from the central s tar (IHinkle. Hall, fc Ridgwayl Il982 



Hinkle. Lebzelter. fc Scharlach 



1997 



Alvarez et al.l 2000 ) a nd complex convective motions i n 



the outer envelope (jPorter. Anderson, fc WoodwardI 119971 : iFreytag. Steffen. fc DorchI |2002| ) . 
These stars have significant mass-loss rates and are often obscured in visible optical bands. 
SiO masers at 43 GHz are however ubiquitous in the extended atmospheres of these ob- 
jects. As compact high-brightness components that are significantly polarized, they act as 
important probes of the astrophysics in the close circumstellar environment, including the 
morphology and relative dynamical influence of magnetic fields. Spectral-line polarization 
VLBI techniques allow direct imaging of this region in Stokes (/, Q, U, V) at a spatial reso- 
lution of /ias, unmatched at other wave-bands. The power of this technique is significantly 
enhanced if synoptic monitoring is conducted over a range of stellar pulsation phase, so that 
the key dynamical drivers can be studied in this region. 

There remain significant uncertainties in integrated astrophysi cal models of the NCSE. 
Stella r pulsation hydrodynamical models are either semi-analvti c (IBertshinger fc Chevalier 



19851 ). or confined to spherical or axi-symmetric numerical studies (lBowerull988l : lBessell. Scholz fc Wood 
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19961 : iHumphreys et al.ll2002l ). Theoretical studies of magnetic fields in AG B stars have pro 



Blackman et al. 



duced a range of predicted magnetic field magnitudes and morphologies 
( 120011 ) have proposed a conventional auj dynamo driven by an inner, differentially-rotating 
stellar core. Photospheric fields of the order of ~ 400 G are predicted in this model, but 
for the case of an isolated (non-binary) AGB star, this model would need a dynamo re- 
supply mechanism to sustain a magn etic field for long enough to shape global mass loss 
(INordhaus. Blackman. fc Frank! 120071). The role of rnagnet ic fields in shaping the global 
mass-loss process is considered bv (IGarcia-Segura et al.lll999l) . In contrast, a convective a^u 
dynamo has been proposed by ISoker k. Zoab J j2002h : in this model, turbulent dynamo fields 
may produce magnetic cool spots 10 — 100 G ) that regulate local dust formation, but have 
no global shaping role in mass- loss JSokei]l2002h . Broader arguments agains t the proposi t 



ion 



4^. ^ 

that global magnetic fiel ds can shape p lanetary nebulae are summarized by lSokerl (120061 ). In 
a numerical MHD study, iDorchl (12004 ) have examined a convective dynamo in a supergiant 
atmosphere and report predicted local fields of up to ~ 500 G. These fields would have 
significant dynamical influence, albeit localized with low filling factors. 

Recent water maser polarimetry of late-type, evolved stars has found evidence for > 10^ 
mG magnet ic fields at the shell radius of the H2Q masers; t his extrapolates to a surface field 
of ~ 10^ G (IVlemmings. van Langevelde. &: Diamond! l2005l ). Recent optical spectroscopy of 
unobscured central stars in planeta ry nebulae (PNe) has measured p hotospheric magnetic 
fields of several kG in these stars (jJordan. Werner, fc O'Toold 12005! ). Magnetic fields of 
this magnitude would be dynamically significant a nd of great astrophy sical importance to 
related studies of collimation mechanisms for PPN (IMeixner et al. !l2004! ). A small subset of 
stellar water maser sources that show linear jets also appear to be magnetically collimated 
( IVlemmings. Diamond, fc Imai!l2006l ). 



As part of a synoptic monitoring campaign of the J = 1 — SiO maser emission toward 
the Mira variable TX Cam at fine time-sampling, this paper presents linear polarization 
images for the first 43 epochs observed, over a pulsation phase range (f) = 0.68 to = 1.82. 
The reduction and analysis of the corresponding total intensity data was published earlier by 
Diamond &: Kemballl (120031 . hereafter Paper I). This paper serves to present the time-series of 
linear polarization images, and provides an analysis of the linear polarization morphology and 
early results on associated proper motions. The associated circular polarization data will be 
presented in a future paper. TX Ca m is a an M8-M10 Mira variable with a pulsation period 
of 557.4 d avs (IKholopov et al.lll985l). and has a Mira period-luminosity distance estimate of 
~ 390 pc (lOlivier. Whitelock. &: Marang||2001l) . It was first observed in total intensity VLBI 



observations by [Diamond et al.! ( 

of the f = {l,2}, J=l — transitions is reported by lYi et al 



I994J ). Recent simult aneous to t al int ensity VLBI imaging 

J2005h . 
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The first polarimetric VLBI imaging of the v = 1, J = 1 — SiO maser emission to 



ward TX Cam was conducted at a single epoch by iKemball &: Diamond! (119971 ). who found 
a predominant orientation of the electric vector position angle (EVPA) tangent to the arc of 
the projected maser shell. This will be referred to as tangential polarization in what follows. 
Subsequent SiO VLBI ir naging polarirn e try o f TX Cam and IRC+10011 at a later single 
epoch was carried out by lDesmurs et al.l (120001 ). who found a similar tangential linear polar- 
ization morphology. A combined optical interferometric and VLBI polarimetric monitoring 
campaign for a sample of 0-rich AGB stars, b ut not includin g TX Cam, is currently under- 



way and is described in a ser i es of papers by ICotton et al. I (120041 ) , ICotton et al. I (120061 ) 



and ICotton. Perrin. &: Lopea (120081 ). Recent single-dish polarimetry of the SiO J = 2 — 1 



transi tion at 86 GHz for a sample of late-type, evolved stars is reported by iHerpin et al. 
J2006h . 



The structure of this paper is as follows. Section 2 describes the observations and data 
reduction methods. The results and discussion are presented in Section 3, and conclusions 
in Section 4. 



2. Observations and Data Reduction 

In this paper we present the polarization reduction of the time-series of 43 synoptic ob- 
servations of the f = l, J=l — SiO maser emission toward TX Cam, previously analyzed 
and reported in total intensity in Paper I. The list of epoch codes, observing dates, and 
associated TX Cam optical pulsation phases are repeated here for reference as Table [H The 
observations in this paper cover the time interval from 24 May 1997 to 19 February 1999, 
corresponding to a TX Cam pulsation phase range 6 = 0.68 to = 1.82. The data were 
observed using the Very Long Baseline Array (VLBAl^) operating in the 43 GHz band, and 
augmented at most epochs by a sing le antenna from the Very Lar^e Arrafl Both tetecope. 
are operated by the National Radio Astronomy Observatory (NRAO§). Auxiliary VLA in- 
terferometer observations were scheduled throughout the period of the VLBA observations, 
to allow absolute calibration of the EVPA of the linearly-polarized SiO emission. Both the 
VLBA and VLA observations and data reduction are discussed in further detail below. Op- 
tical variability data, showing the pulsation phases in relation to the scheduling of these 

^http:/ /vlba. aoc.nrao.edu 
^http:/ /via. aoc.nrao.edu 

•^The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated 
under cooperative agreement by Associated Universities, Inc. 
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VLBA observations, is shown in Figure [H reproduced here from Paper I. 

2.1. VLBA Observations and Reduction 

The data were recorded at each antenna in dual-circular polarization in two 4 MHz 
baseband spectral windows, each digitally sampled at the full Nyquist rate of 8 Mbps in 
1-bit quantization. The lower spectral window was centered at a fixed topocentric frequency 
corresponding to the v = 1, J = 1 — SiO transition, at an assumed rest frequency uq = 
43.122027 GHz and a systemic velocity Vlsr = +9 km s~^toward TX Cam, Doppler-shifted 
to the center of the array and the midpoint of the observations at each epoch. No further 
real-time Doppler corrections were applied; these were performed in post-processing. For the 
majority of the VLBA epochs in Table [1], the v = 2, J = 1 — transition was simultaneously 
observed, and centered in the second spectral window; the analysis of these data will be 
presented in a future paper 

Each VLBA observing epoch was scheduled in a comparable local sidereal time (LST) 
range, as permitted by telescope operations. Three fixed schedule formats were used over 
the course of these observations. For epochs {A-0} a schedule of 6-hour duration was used, 
balanced between the target source, TX Cam (3.5 h), the primary continuum calibrator, 
J0359+509 (1.1 h), and secondary calibrators 3C454.3, J0609-157, and 3C286 (each 0.22 
h). A single scan on 3C286 was included initially to explore an alternative absolute EVPA 
calibration method, as discussed further below. This was not successful however, and for 
epochs {P-X} the 3C386 scan was re-allocated to TX Cam (3.6 h), and the schedule re- 
balanced. For epochs {Z-AQ}, an 8-hour schedule format was used, but 2 hours of each run 
were shared with a 22 GHz project, leaving the total 43 GHz observing time, and its division 
between sources, unchanged. 

The data were correlated at the VLBA correlator in Socorro, NM, adopting a field 
center position for TX Cam of (ajaooo = 05^00™51M86, Sj2ooo = 56'^10'^54^341). The 
correlator accumulation interval At^. was constrained by output data rate limits in place at 
the correlator at that time. A value Ate = 6.03 s was used for epochs {A-G,LK}, and reduced 
to Ate = 4.98 s for epochs {H,M-AR}. All polarization correlation products (RR,RL,LR,LL) 
were formed in each correlation accumulation interval from the dual-circular data acquired 
at each antenna. In this spectral-line polarization mode the VLBA correlator has a channel 
count limit of Nchan = 128. This produced auto- and cross-power spectra in each 4 MHz 
baseband with a nominal velocity spacing of ~ 0.2 km s~^. 

The measured visibility data were reduced using spectral-hne polarization calibration 
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methods described by lKemball. Diamond, fc CottonI (Il995l ) and lKemball fc Diamond ( 19971 ). 
as implemented in a semi-automated pipeline using a modified version of the AlPq^ pack- 
age maintained by the NRAO. These analysis methods are optimized for the reduction of 
spectral-line polarization VLBI observations of the type described here; in particular they 
make no implicit assumption that Stokes = 0. The SiO molecule is non-paramagnetic, 
and the mean degree of circu lar polarization in the SiO maser rotational trans itions is small 
but not zero, rric ~ 1 — 3% (IBarvainis et al.l 119871 : iKemball fc Diamond! 119971 ) . In the data 
reduction method employed here, all antenna-based group delay, fringe-rate, and phase cali- 
bration, including final phase self-calibration, was performed relative to a reference antenna 
in a reference receptor polarization. This calibration was transferred to the orthogonal re- 
ceptor polarization by solving for all relevant phase and group delay offsets using continuum 
calibrator observations. This generalized reduction method retains positional coincidence 
between Stokes / and V images, while making no implicit assumption that Stokes V is 
identically zero (IKemball. Diamond, fc Cottonlll995[ ). Amplitude calibration was performed 
with the same goal of preserving the low Stokes V signatur e. Each circ u lar re ceptor po- 
larization was calibrated independently using the method of iReid et al. I ( 1l980l ) to fit the 
autocorrelation spectra in each receptor polarization to a template reference spectrum in 
that polarization. The template spectra in each receptor polarization were calibrated on 
an absolute scale using system temperature measurements reported during the observations 
and known point-source sensitivities and gain curves published for the VLBA antennas. A 
final single differential polarization gain correction was determined by cross-fitting the spec- 
tra between receptor polarization at the template scan to measure a relative scale factor. 
For the low integrated rric known for SiO masers, this approximation is warranted for the 
integrated autocorrelation spectrum; it would not be appropriate for individual cross-power 
spect ra, for example, as some S iO maser components can be substantially circularly polar- 
ized (IKemball fc Diamond! 119971 ). This final differential polarization gain correction adjusts 
for errors in the reported system temperatures and point source sensitivities, as used in the 
a priori absolute calibration of the template spectrum in each receptor polarization. We 
estimate these uncertainties very approximately as 5-10% for the VLBA at this frequency. 
We do not correct for atmospheric opacity at the pointing position and time of observation 
of the template spectrum for the data reported here. However, this has no polarization 
dependence, and accordingly does not affect our measurement of linear or circular polariza- 
tion percentage or orientation; it only introduces a modest uncertainty in our absolute fiux 
density scale. This uncertainty is mitigated by the fact that the template spectrum is chosen 
for its signal-to-noise ratio (SNR) and peak magnitude; these factors ensure that it is almost 



"'http: / /www. aips.nrao.edu 
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always at a sufficiently high elevation to minimize atmospheric opacity contributions. 

As phase self-calibration is used, our final VLBA images of the SiO maser emission 
toward TX Cam retain no absolute astrometric information or alignment. In Paper I, ze- 
roth frequency-moment total intensity images from successive epochs were initially aligned 
using spatial cross-correlation. This initial alignment was then refined using feature-based 
registration in which the deviation in proper motion paths predicted for individual maser 
components was minimized. For the polarization and total intensity images produced and 
reported in the current reduction, we adopt the alignment frame developed in Paper I as 
an absolute reference. Each of the zeroth-moment total intensity images produced in the 
current work was registered relative to its corresponding epoch in the Paper I frame using 
Fourier-based matched filtering. A robust, weighted-sum estimator of the filter output peak 
was used to measure the two-dimensional spatial offsets. This was found to out-perform an 
elliptical Gaussian fit to the peak, likely because it is less sensitive to small morphological 
differences in the total intensity images being aligned. Based on manual cross-checks of the 
positions of individual SiO maser components, we believe the registration of our polarization 
images relative to the alignment frame determined in Paper I has an accuracy of lOOfias or 
less. Note that the images in Paper I used to produce the Stokes / movie published there 
were enhanced in contrast by applying the transformation a/7 and a subsequent threshold 
cutoff vT < 0.7. For consistency, the same transformation was applied to the total intensity 
images reduced here before the matched filter alignment. 

The calibrators, J0359+509, 3C454.3, and J0609-157 were reduced separately using 
co ntinnuum polarization VLBI techniqu e s for millimeter wavelength observations described 



by iKemball. Diamond, fc Pauliny-TothI (119961 ). As the only calibrator source with suffi 



cient parallactic angle coverage, continuum instrumental polarization leakage terms were 
determined from the observations of J0359+509 alone and applied to the other calibrator 
sources. Continuum polarization images in Stokes (/, Q, U) were then produced for each 
calibrator at each epoch. As noted above, the calibrators 3C454.3 and J0609-157 were only 
observed for a single 13-min scan at each epoch, primarily as visibility calibrators, and 
have limited image sensitivity and poor wv-coverage as a result. The image quality for these 
sources was accordingly often poor, especia lly for 3C454.3, which has more complex structure 



( IKemball. Diamond. &: Pauliny-Tothlll996l ) than the more compact J0609-157 



2.2. VLA Observations and Reduction 



The absolute EVPA of any linearly-polarized emission measured in our VLBA observa- 
tions is unknown as no instrumental measurement of the absolute R-L phase difference at the 



- 8 - 



reference antenna is possible. Although the VLA similarly lacks the capability to make this 
measurement, it can provide absolute astronomical calibration of EVPA relative to a small 
subset of primary polarization calibrators given the lower spatial resolution of the array. 

We have adopted the traditional approach of calibrating the EVPA of our VLBA ob- 
servations by using the VLA to transfer absolute polarization calibration to these data via a 
sufficiently compact secondary polarization calibrator suitable for imaging with both arrays. 
We use the continuum source J0359+509 as this transfer calibrator. Short (1-2 h) VLA 
observations were accordingly requested across the time-span of our VLBA observations, 
and were scheduled when logistically possible on the VLA. The dates of the auxiliary VLA 
observation associated with our TX Cam monitoring campaign are listed in Table [2j 

In these VLA observations we measured the absolute EVPA of J0359-f-509 relative to 
the primary polarization calibrator 3C138, which is the only primary polarization calibrator 
available in this right ascension range. Our observations were primarily scheduled in contin- 
uum mode, observing with two 50 MHz spectral windows and full polarization. Instrumental 
polarization calibrators were chosen separately for each VLA run, as it was not possible to 
schedule these observations in identical local sidereal time ranges. Each observation included 
J0359+509 and 3C138 however, in order to address the primary scientific goal. 

Two dominant sources of systematic error affected these observations: i) their short 
duration, which at times limited the fidelity of the instrumental polarization calibration due 
to the reduced parallactic angle range; and ii) the extended structure of 3C138 at this angular 
resolution and observing frequency, making it a sub-optimal polarization reference in the 
long-baseline VLA configurations. The deficiencies in 3C138 notwithstanding, it remained 
our only viable primary polarization calibrator in this LST range. 

The data were reduced using standard heu ristics for Q-band V LA observations, as de- 



scribed in an appendix of the AIPS cookbook (iGreisen et al. 1120071 ). During data loading, 
an opacity correction was applied using a weighted estimate of the zenith opacity based 
on the surface weather, with a relative weight of 0.75, and a seasonal model, weighted by 
0.25. A standard gain curve supplied for the VLA was used. Post-facto antenna-position 
corrections were then applied as needed. All subsequent reduction was performed for each 
50 MHz spectral window independently. The data were edited interactively and the primary 
fiux scale established for 3C138 using the model provided within the AIPS package. An 
initial phase calibration over a solution interval of 20 s was performed, using an a priori 
model for 3C138 (to mitigate source resolution) and any uv-range limits recommended for 
other calibrators in the schedule, but which were otherwise assumed to be point sources. An 
amplitude gain correction was then determined for each antenna, scan, and receptor polar- 
ization on each continuum source. This allowed the fiux density of the continuum calibrators 
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other than 3C138 to be determined relative to the prior flux density scale established for 
3C138. A joint solution for instrumental polarization was performed across all continuum 
calibrators, and final images in Stokes (/, Q, U) produced for each continuum source in each 
50 MHz spectral window separately. An image pixel spacing of ^th of the expected spatial 
resolution in each VLA configuration was adopted, assuming that the VLA configurations 
have the following spatial resolution s at 43 GHz: A) 0.05"; B) 0.15"; C) 0.47"; and, D) 



1.5" (jUlvestad. Perley. fc Tayloiil2007l ). For hybrid configurations, such as AB, the geometric 
mean of the nominal resolution in each configuration was used. The images were cleaned 
to the la thermal noise limit predicted for each source based on the number of calibrated 
visibilities, time and frequency averaging, and known system equivalent flux density. 

The absolute EVPA of J0359+509 was determined separately in each 50 MHz spectral 
window, relative to the value measured for 3C138 in that window, here assumed to have 



an absolute EVPA xevpa = -14° (IPerlev fc Tavloil 120031 ). The final absolute EVPA x of 



J0359+509 was computed as the mean of the EVPA measured independently in each 50 
MHz spectral window, here denoted by subscripts 1 and 2, with an estimated standard error 
= ■^\xi~X2\- This error estimate is approximate due to the small number of independent 
measurements. The measured absolute EVPA for J0359+509 is listed in Table [2] and plotted 
in the top panel of Figure [2l A weighted cubic polynomial fit over a scrolling window five 
samples wide is drawn as a solid line in the upper panel. This fit was used to interpolate the 
absolute EVPA measured with the VLA to the dates of the VLBA observations, as listed 
in Table [TJ We note however that the VLA data are sparsely-sampled in time, especially 
during a sharp rise in EVPA in late 1998. Continuum extra-galactic sources are significantly 
variable in linear polarization, and this interpolation scheme does not capture all source 
variability at the current sampling. 

The EVPA measured for J0359+509 in the VLBA observations is plotted in the mid- 
dle panel in this figure. This measured EVPA is not absolutely calibrated and is in- 
stead uncertain by an additive rotation equal to half the RCP-LCP phase difference at 
the VLBA reference antenna, which in the VLBA reduction was chosen to be Los Alamos. 
Error bars cr^were estimated approximately for the VLBA continuum polarization images 
from the measured off-source rms, ag and au, in Stokes Q and U respectively as = 

2{Q'2+u'^ \/'^uQ'^ + ^V^'^ \J^^ where is the number of pixels over which (Q, U) were mea- 
sured, and aq and au were corrected for noise correlation resulting from c onvolution wit h 



the restoring beam (of mean FWHM ~ 300;Ltas) using the relation given by ICondonI (jl99ll ). 



The predicted EVPA rotation for each epoch, obtained as the difference between the 
interpolated VLA and measured VLBA EVPA at each epoch, is plotted as a series of points 
in the lower panel. We expect, from instrumental considerations, that the RCP-LCP phase 
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difference at a VLBA antenna should be relatively stable over long periods, barring any 
changes in receivers or electronics. We note the bimodal distribution of points that has a 
vertical break coinciding with the only receiver change at the reference antenna (Los Alamos) 
during the period of these observations. This occurred on 30 April 1998. Accordingly, we 
fit a cubic polynomial to the data separately on either side of this break to derive the final 
EVPA rotation per epoch. This was the lowest-order polynomial to yield a reasonable fit 
to the data. The significant deviations in late 1998 are likely due to the systematic under- 
sampling of the VLA EVPA measurements during this period, as noted above. We estimate 
the peak-to-peak error in absolute EVPA alignment at ~ 10° — 20° using this method. 



3. Results and Discussion 

The total and linearly-polarized intensity images of the v = 1, J = 1 — SiO emission 
toward TX Cam over the pulsation phase range covered by this paper are shown in Figures [3] 
through O In these composite figures, the enhanced-contrast Stokes I image is plotted as 
a background contour plot, overlaid by vectors representing the linearly-polarized emission. 
The vectors are drawn with a length proportional to the linearly-polarized intensity and 
at the position angle (N through E) corresponding to the EVPA of the linearly-polarized 
emission. These figures are the zeroth-moment (over frequency) of the associated Stokes 
(/, Q, U) brightness in the parent image cubes at each epoch. 

These polarization images are presented in an animated sequence in the electronic mate- 
rials associated with this paper, in the form of a movie of the linear polarization morphology 
evolution of the f = l, J=l — SiO maser emission toward TX Cam over the pulsation 
phase range = 0.68 to = 1.82. 

The phase range covered is the first pulsation cycle of our longer synoptic monitoring 
campaign. Several striking properties are evident from this time-series of figures. Over 
this pulsation phase range, the linear polarization EVPA is predominantly tangential to the 
projected total intensity SiO maser shell. In addition, for the majority of the SiO maser 
components, the EVPA of individual features appears to persist over their component life- 
times. Finally, the linearly-polarized intensity of the maser emission is generally brightest 
at the inner projected shell boundary and decreases in magnitude at larger radii. 

The dominant tangential li near-polarization morphology is consistent with that reported 



by iKemball &: Diamondl (119971 ) at a single epoch, in the first VLBI imaging polarimetry of 



this SiO transition toward TX Cam. The earlier observations were undertaken on 1 Dec 1994, 



so precede the data in the current paper by two pulsation cycles. iDesmurs et al.l (120001 ) 
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confirmed the tangential linear polarization morphology toward TX Cam in observations 
in April 1996. Taken together with the current data, this suggests that this morphology 
may be an inter-cycle property of the SiO maser emission toward this star. It is not yet 
clear however, whether this is a generic pr operty of SiO li n ear p olarization in this transition 
toward late-type evolved stars in general. iDesmurs et al.l (120001 ) report a similar tangential 
morphology in the SiO maser emission toward IRC-l-10011. However, for a sample of Mira 
variables monitored inf = l,2, J=l — SiO maser emission in a full-polariz a tion V LBI 
imaging study, and reported in a series of papers by ICotton. Perrin. fc LopezI (|2008| . and 
references therein), t he results are more m ixed. Early images suggested no pervasive linear 
polarization pattern (jCotton et al. II2004I ): however, later, a bimodal EVPA distribution (in 
the sense of bei ng either parallel or perpendicular to the projected shell) was reported for U 



Ori and o Get (ICotton et al. 1120061) . However, not all sources for which images are shown 



by ICotton. Perrin. fc LopezI (120081 ) show tangential linear polarization structure. Further 
observations of larger samples are needed to resolve this important question. 

There are several possible theoretical explanations for tangential linear polarization 
structure in SiO maser emission, as we see in the data presented here. However, these inter- 
pretations are conditional on the theoretical assumptions made concerning the underlying 
maser emission; we focus here only on the observational implications of these theoretical 
issues in what follows. 

Silicon monoxide is a non-paramagnetic molecule and t he rotat i onal t ransition consid- 
ered here is in the small-splitting Zeeman regime. Foll owing Elitzurl (1996 ), the ra tio of 
the Zeeman splitting to the Doppler linewidth Avd is (IKemball fc Diamondl 119971 ): 



xb = 8.2 X 10" 



km.s~" 



-1 



For the case of <C 1, as equation [T] indicates is applicable t o SiO, compe t ing models for 



the transport of polarized maser emission have been presented (lElitzurlll996l : IWatsonI 12002 



and references therein). Both seri es of papers extend the paramete r space and generalize 
the original foundational work of iGoldreich. Keeley. fc KwanI (119731 ). There also remains 
uncertainty in the literature as to whether SiO masers are pumpe d by collisional (lElitzur 
1980|) or radiative mechanisms (iBujarrabal fc Nguyen-0-Rieulll98ll ). Finally, an additional 
consideration is the degree to which m— anisotropic pumping of the magnetic substates 
dominates. We neglect for the purposes of this discussion the role of magnetic or velocity 
gradients in the masing regions, and any considerations that apply specifically to transitions 
with spins higher than J = 1 — 0. 



At the outset we need to determine in which theoretical regime of maser polarization 
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propagation our observational data reside. Specifically, there are several key parameters that 
determine this theoretical interpretive framework. These include the relative magnitude of 
the stimulated emissi on rate, R, the cumulative coUisional and radiat ive decay rate, F, and 
the Zeeman rate, gfl ( IGoldreich. Keeley. fc Kwanlll973l : IWatsonll2002l ). 



To estimate R, we adopt equation 1 of iPlambeck. Wright, fc Rad (120031 ) to derive the 
stimulated emission rate of the f = l, J=l — SiO transition as: 



R=23 



2 X 1010 K 



dQ 



10- 



sr 



(2) 



where Tg is the maser brightness temperature and dQ is the estimated maser beaming 
angle. The relative magnitude of R and F defines the degree of maser saturation, with 
saturation increasing for increasing ^. The vibrational radia tive transitions v (v — 1) for 
low- J rotation levels have radiative decay rates F ~ 5f (IKwan &: Scovillelll974j : lElitzur 



19921 ). These radiative decay rates dor ainate o v er SiO -H2 or SiO-H coUisional de-excitation. 
From the coUisional rate equations in lElitzurl (jl980l ). we calculate F < Is"^ for collisions, 
adopting a representative temperature T=1300 K and atomic hydroge n density uh = 5 x 10^° 



cm ^ in the Si O maser region, as inferred fr o mlReid fc MentenI (119971). The Z eeman rate for 



SiO is given by iNedoluha fc WatsonI (Il990f ): IPlambeck. Wright, fc Raol (l2003f ) as: 



2% 



200 



-1 



(3) 



In the limiting regime gQ 3> ^ F, and assuming m-isotropic pumping, the IGoldreich. Keeley. fc Kwe 



( 119731 ) solutions for the linear polarization are applicable. In this case the fractional linear 
polarization depends on the angle between the line of sight and the magnetic field, and 
takes the form g(9) = f = "^'.l^y for 6 > Ob and q{e) = I ioi 9 < Ob, where Ob ~ 35° 
(IGoldreich. Keeley. fc Kwaru 1 19731 : lElitzurl Il992l ). At a break angle. Op ~ 55°, defined as 



the point where sin^ 6*^? = |, and Stokes Q changes sign, the observed EVPA on the plane 
of the sky switches from parallel to the project ed magnetic field {9 < 9f) to pe r pendic- 
ular to the projected magnetic field {9 > 9p) (IGoldreich. Keeley. fc KwanI Il973l : lElitzur 
1992[). We note that th i s cha nge of EVPA over 90° has been detected observationally by 
Vlemmings fc Diamond! (120061 ) for a water maser component in the source W43A. 



The degree of saturation ^ required to achieve the the limiting magnitude of fractional 
linear polariz ation differs however between the two th eoretical models of maser polarization 
propagation (jNedoluha k. WatsonI Il990l : lElitzurl Il996l ) . A satur ation of ^ ~ 30 i s requ ired 
to achieve 70% linear polarization in the models summarized by lWatson fc Wyldl (120011 ). In 



- 13 - 



the model of lElitzurl (jl996l ). however, the full polarization solution is realized well before 
saturation. 

In addition, the observational interpretation of linearly-polarized SiO maser emission 
is further influenced by whether the maser inversion is assumed to be due to collisional 
or radiative pumping. If the pumping is assumed to be radiative, then the effects of 
m-anisotropic pur nping of the magnetic sublevels m of the rotational transitions needs 
to be considered (IBuiarrabal fc Nguven-Q-Rieul Il98ll ) along with related radiative effects 
(iRamos. Degl'Innocenti. fc Buend l2005l ). In the case of stellar SiO masers, there is an 
inherent anisotropy in the angular distribution of the photons from the central star inci- 
dent on the masing region, and, in this model, m-anisotropic pumping can produce large 
l inear polarization magni t udes, oriented tangen tially to the incident radial stellar photons 
(jWestern fc Watson! Il983l : iDesmurs et al.l 120001 ) . For m-anisotropic pumping, the observed 
EVPA is still likely to be either parallel or perpendicular to the projected m agnetic field o n 
the plane of the sky, if the conditions gQ > R and gQ > T hold, even weakly (jWatsonll2002l ). 
however the shape of the function q{6) and the break angle 6f are likely to differ from the 
standard Zeeman case above. If the magnetic field is sufficiently weak so that gQ ~ R holds, 
then the E VPA need no l onger be strictly parallel or perpendicular to the magnetic field in 
this model (IWatsonll2002l ). Given the degree of order we observe in the linear polarization 
morphology, the data suggest that gfl ^ R in the current study. 

Both theoretical models of maser polariz ation predict circular polarization tha t is broadly 
proportional to the magnetic field strength (jWatson &: Wyldll200ll : lElitzurl Il996f ). but with 
differences in predicted line-shape and dependence on saturation as well as the angle 6 
between the magnetic field and the line of sight. Given xb for SiO ma sers, circular polariza- 
tion at the several percent level predicts Gauss-level magne tic fields f Barvainis et al.lll987l : 
Kemball fc Diamondl 119971 ). In the models summarized by IWatsonI (120021 ). circular polar- 
ization can also arise from non-Zeeman effects caused by a difference between the optical 

axes along the maser propagation path and the direction of linear polarizatio n. Circular 

polarization can be produced in this model when R f» gfl (iNedoluha fc Watsonlll990l ) or by 
turbulent rotation of th e magnetic field along the line of sight in the regime gfl ^ {R, T} 
Jwiebe fc Watson! Il998h . 



To constrain the theoretical regime in which our current data reside, it is first useful to 
estimate the relative magnitude of gfl, R, and F. There are two main sources of uncertainty 
in this estimation, namely the strength of the magnetic field, B, and the unknown maser 
beaming angle dfl used to calculate the stimulated emission rate in equation [2] (equivalently, 
the degree of saturation). We make an initial estimate of B by comparing energy densities. 
A magnetic field B implies a magnetic energy density: 
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2/zo 



4 X 10" 



(4) 



For our adopted uh = 5 x 10^° cm~^ and T = 1300 K, the thermal energy density is 
^nnkT ~ 1.3 x 10^'^ J.m^^, and the bulk kinetic energy density is ^pv"^ ~ 2.0 x 10^^ J.m~'^, 
where p is density and we adopt a bulk motion v = 7 km s~^from Paper I. Equating these 
energy densities suggests a zeroth-order magnetic field estimate in the range ~ 580 — 720 
mG. For B ~ 0.6 G, equation [3] predicts |^ ~ 120 s~^. Tangential shock compression driven 
by stellar pulsation could increase th e magnetic field strength by an order of magnitude 
JSokei]l2002l : iHartquist fc Dvsonlll997[ ). Adopting ~ 23 s-\ F ~ 5 s"\ and gQ > 10^-^ 
would place us in the regime: gfl ^ R ^ T or gQ > R ^ T. In this regime, and absent 
significant Faraday rotation, the standard Zeeman interpretation of the parallel or perpen- 
dicular orientation of the EVPA relative to the projected magnetic field would apply, and 
the linearly-polarized morphology could either be tracing a globally-organized stellar mag- 
netic field, or a sparse strong field that is only dynamically significant on local scales. There 
are arguments against a globally organized stellar magnetic field in isolated (non-binary) 
TP-AGB stars, primarily the likely low rotation velocity at this point in their evolutionary 
path in the H-R diagram, and the difficulties of sustaining a regular a u dynamo in stars 
with such an extended mantle, amongst other considerations (ISokerll2006l ) .However, as noted 
above, there are several strands of o bservational evidence supporting surface magne tic fields 
of order 10^"^ G in Mira variables ([Vlemmings. van Langevelde. fc Diamondl 120051 ) and the 
central stars of PNe (jjordan. Werner, fc O'Tooy 2005h respectively, as well as rn agnetic jet 
collimation in post-AGB transition objects (j Vlemmings. Diamond, fc Imail 120061 ). 



A sparse field in an AGB extended atmosphere could arise from an a'^uj convective dy- 
namo, as noted above, and might incl ude local loops rising above the photosphere, ot her flar- 
ing events, and magnetic cool spots (ISoker fc Claytonlll999l : ISoker fc Kastnerl l2003l ) . These 
localized fields are also expected to be shock- compressed tangentially due to the outwardly- 
propagating shock waves from the pulsating central star, introducing a preferred geometry 
for field enhancement. This is a complex three-dimensional picture, but it is also important 
to remember that SiO masers sample only those regions of the NCSE that lie relatively 
close to the orthogonal plane to the line of sight. This constraint arises from the observed 
tangential amplification of the total intensity emission. 

If however, the masers are pumped radiatively and the radiative pumping is m- anisotropic, 
then the tangential linea r polarization morphology could origir iate from that mechanism 
(IWestern fc Watsorull984j ) , even in the non- magnetic case gfl ~ (IRamos. DegFInnocenti. fc Bueno 
2OO5I ). For the magnetic case, gil > 0, the Hanle effect can, in certain cases, produce an 
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EVPA rotation exceeding 45 deg relative to the tangential direction (IRamos. Degl'Innocenti. &: Bueno 



20051 ). In the case of m- anisotropic pumping in general, we would also expect to measure a 



functional shape for q(9 ), and associated break a ngle 6^, that differs from the standard Zee- 



mag nitude JWatsonI 12002 ). 



man form n oted above (IWestern fc Watsonlll984j ). and a higher degree of linear polarization 



In Figure [TO] we examine the variation of mean fractional linear polarization magnitude 
frii{(j)) against pulsation phase. This value was computed at each epoch from the sum of 
Stokes / and P = ^/Q^TTP, using a pixel mask defined as the set of all pixels with a 
total intensity exceeding three times the Stokes / image rms, and corrected for statistical 
bias in P. The value of riii{(f)) has a minimum near ~ 1 and a maximum near ~ 1.5. 
It has higher values when the inner edge of the shell is well-defined and dominates the 
overall total intensity morphology. Unfortunately this result can be interpreted either as 
evidence for linear polarization driven primarily by magnetic field enhancement through 
shock compression, or is supportive of m-anisotropic pumping as the dominant origin of 
linear polarization: both are expected to be stronger in the inner shell boundary. We note 
that the mean fractional linear polarization is broadly inversely correlated with the shell 
width. 

We identify several key uncertainties that need to be further constrained by observa- 
tions before the remaining theoretical polarization transport uncertainties can be resolved, 
and that are amenable to observational tests: i) the statistical incidence of tangential linear 
polarization morphology in Mira variables as a class; ii) independent measurement of the 
exact degree of saturation of stellar SiO masers; and iii) detailed measurements of the func- 
tional form of q{6) in maser regions where this is possible. We plan to address these issues 
in future papers. As yet, there is insufficient observational data to distinguish the relative 
influence of m-anisotropic pumping on the magnitude or orientation of the linearly-polarized 
SiO emission 

The linear polarization maps we present here have signiflcant fine-scale, locally-ordered, 
polarization structur e , as n oted in the first SiO polarization observation s of TX Cam by 



Kemball fc Diamondl (119971 ). Recently, in their sample of 0-rich AGB stars. ICotton. Perrin. fc Lopez 
(120081 ) have cited locally-ordered polarization morphologies along jet-like features pointing 
radially toward the central star. Due to the persistence of the SiO linear polarization prop- 
erties over a sufficient fraction of the component lifetime in the polarization movie presented 
here, we are able to explore the time evolution of the polarization morphology for individual 
maser components. In this paper we undertake a preliminary analysis confined to two iso- 
lated components selected subject to the criteria that they: i) follow continuous trajectories 
across a sufficient range of epochs; ii) are located at a larger projected separation from the 
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inner shell boundary, so as to minimize the possible impact of strong anisotropic pumping 
effects; and, iii) have compact structure and high SNR. 

In Figure [TT] we show this analysis for a prominent infalling component at the cen- 
ter north of the projected shell, that is visible over epochs {S-AH}. In Figure [TTl the 
linear polarization structure of this northern component is shown at the twelve epochs 
{T,U,V,W,X,Z,AA,AB,AC,AD,AE,AG} spaced along the trajectory over time, using the 
same composite representation of linearly-polarized emission as used in Figures [3] to [HI but 
with a more finely-sampled total intensity contour scale, described in further detail in the 
figure caption. The outermost epochs {S,AH} of the component trajectory are omitted due 
to the low SNR when the component appears and disappears. In common with Figures [3] to 
[HI each panel in Figure [TT] is a zeroth-moment image over frequency; the pronounced com- 
ponent elongation visible in these zeroth-moment component images at several epochs along 
the trajectory is most readily explained in terms of a velocity gradient along the elongation 
major axis. These velocity gradients are most directly assumed to arise from linear accel- 
eration along a three-dimensional component trajectory inclined to the line of sight to the 
observer. The weighted-mean component positions, with associated error bars, are plotted 
in Figure [T^ for each epoch along the trajectory. At each component position in this figure, 
a bold vector of uniform length (over epoch) is drawn at the orientation of the integrated 
component EVPA. A least-squares fit to the comp onent positions , incor porating x- and y- 



errors, was performed using software described by IWeiner et al. I (120061 ). The fitted mean 
projected trajectory is extrapolated in Figure [H] across the full diameter of the SiO maser 
shell at epoch {Z}, here chosen as an epoch near the mid-point of the northern component 
trajectory shown in Figure [12l Figure [H] shows that the projected component trajectory is 
non-radial with respect to the likely stellar position relative to the SiO shell distribution; 
this qualitative assertion remains valid within the statistical errors of the projected linear 
fit parameters. If dynamically-significant gas pressure gradients are not present, then purely 
gravitational infall would be radial toward the central star. We note that this component 
has an EVPA oriented predominantly orthogonal to the projected proper motion vector. 
As noted above, the location of this component is at a significantly larger radius than the 
inner shell boundary, and therefore minimizes the possibility of strong anisotropic pumping 
effects. Over parts of the earlier epoch range {C-N}, a component is visible approximately 
3.25 mas to the east of the northern component considered here. This north-eastern com- 
ponent is at a significantly lower SNR, and is not clearly detected in linear polarization at 
epochs {G-I,M-N} as a result. Given the proximinity to the northern component, however, 
we plot the counterpart of Figure WI\ for this component as Figure [T31 This shows the north- 
eastern component moving outward over these epochs with an EVPA again predominantly 
orthogonal to the projected proper motion vector. 
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A second isolated component is visible over epochs {AA-AL} to the east of the mid- 
latitude of the projected shell, and we plot the counterparts of Figures [TTIfT2] for this eastern 
component as Figures [TMT^ This component is infalling over this epoch range and the 
EVPA is primarily either parallel or perpendicular to the projected proper motion vector. 

This preliminary polarization proper motion analysis for the initial components consid- 
ered here, which were chosen for their location in the outer shell where radiative anisotropic 
pumping effects are less likely to predominate, show two key properties: i) the inferred 
magnetic field orientation is likely to be parallel to the proper motion vector, as the EVPA 
orientation is either parallel or perpendicular to this proper motion; and ii) the motion is 
sometimes non-radial relative to the central star, in contrast to what is expected for ballistic 
motion in the stellar gravitational field. These properties are consistent with component 
motion along magnetic field lines in the extended atmosphere. An alternative explanation 
is that a local magnetic field is being dragged and elongated by the gas motion, and the 
component motion is not governed dynamically by the magnetic field. This explanation is 
not excluded however, but we believe it is not favored by the preponderance of the cur- 
rent observational evidence. This issue is also discussed in connection w ith jet-like feature 
extensions in the SiO images obtained by (jCotton. Perrin. fc Lopeal2008l ). 



We also observe several 90-degree EVPA reversal at the inner-shell b oundary, as noted 
in the first TX Cam polarization imaging flKemball &: Diamond! 119971 ). These features 
are particularly prominent in the south-west edge of the shell, from epochs {AJ-AR}. 
These EVPA reversals can b e exp lained either as transitions across the break angle in the 
Goldreich. Keeley. fc KwanI (119731 ) so lutions or, in the case of m-anisotrop ic pumping, a 
transition in the degree of anisotropy (IRamos. DeglTnnocenti. fc Buendl2005l ). These issues 
will be considered in more detail in a future analysis. 



4. Conclusions 



We have presented a finely-sampled time series of 42 images of the linearly-polarized 
V = 1, J = 1 — SiO maser emission toward the Mira variable, TX Cam over a pulsation 
phase range (p = 0.68 to = 1.82. These data augment a total intensity aniraated s equence 
previously published for this pulsation phase range by [Diamond fc Kemballl (120031 ). From 
our initial analysis we conclude: 



1. The predominant linear polarization morphology is one in which the EVPA of the 
linearly-polarized SiO emission is predominantly tangential to the projected shell of 
SiO maser emission in total intensity. This morphology is observed for all epochs 
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presented here. 

2. The hnear polarization EVPA of individual SiO maser components is relatively stable 
over their component lifetimes for many individual components visible in the images. 

3. The magnitude of the linear polarization appears brighter at the inner region of the 
shell boundary, and decreases for larger radii. 

4. An analysis is presented that argues that we are in the theoretical regime gQ:$> R^T 
or gD, > R ^ r for maser polarization propagation. There is, as yet, insufficient 
observational data to distinguish the relative influence of m-anisotropic pumping on 
the magnitude or orientation of the linearly-polarized SiO v — 1, J — 1 — emission. 

5. An initial polarization proper motion study is presented for two outer maser compo- 
nents that provides some supporting evidence to the argument that magnetic fields are 
dynamically significant in the extended atmospheres of late-type, evolved stars. 
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Table 1. Observing dates and epochs 



Epoch code Observing date Optical phase ^ 

(0) 



BD46A 


1997 


May 24 


0. 


,68 


± 


0, 


.01 


BD46B 


1997 


June 7 


0. 


,70 


± 


0, 


.01 


BD46C 


1997 


June 22 


0. 


,73 


± 


0, 


.01 


BD46D 


1997 


July 6 


0. 


,75 


± 


0, 


.01 


BD46E 


1997 


July 19 


0. 


,78 


± 


0, 


.01 


BD46F 


1997 


August 2 


0. 


,80 


± 


0, 


.01 


BD46G 


1997 


August 16 


0. 


,83 


± 


0, 


.01 


BD46H 


1997 


August 28 


0. 


,85 


± 


0, 


.01 


BD46I 


1997 


September 12 


0. 


,87 


± 


0, 


.01 


BD46J 


1997 


September 26 


0. 


,90 


± 


0, 


.01 


BD46K 


1997 


October 9 


0. 


,92 


± 


0, 


.01 


BD46L ^ 


1997 


October 25 


0. 


,95 


± 


0, 


.01 


BD46M 


1997 


November 8 


0. 


,98 


± 





.01 


BD46N 


1997 


November 21 


0. 


,00 


± 


0, 


.01 


BD460 


1997 


December 5 


1. 


,03 


± 


0, 


.01 


BD46P 


1997 


December 17 


1. 


,05 


± 


0, 


.01 


BD46Q 


1997 


December 30 


1. 


,07 


± 


0, 


.01 


BD46R 


1998 


January 13 


1. 


,10 


± 


0, 


.01 


BD46S 


1998 


January 25 


1. 


,12 


± 


0, 


.01 


BD46T 


1998 


February 5 


1. 


,14 


± 


0, 


.01 


BD46U 


1998 


February 22 


1. 


,17 


± 


0, 


.01 


BD46V 


1998 


March 5 


1. 


,19 


± 


0, 


.01 


c 


1998 


March 21 


1. 


,22 


± 


0, 


.01 


BD46W 


1998 


April 07 


1. 


,25 


± 


0, 


.01 


BD46X 


1998 


April 19 


1. 


,27 


± 


0, 


.01 


BD46Z 


1998 


May 10 


1. 


,30 


± 


0, 


.01 


BD46AA 


1998 


May 22 


1. 


,33 


± 


0, 


.01 


BD46AB 


1998 


June 6 


1. 


,35 


± 


0, 


.01 


BD46AC 


1998 


June 18 


1. 


,37 


± 


0, 


.01 
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Table 1 — Continued 



Epoch code Observing date Optical phase 

(0) 



BD46AD 


1998 


July 3 


1 


40 


± 





01 


BD46AE 


1998 


July 23 


1 


44 


± 





01 


BD46AF 


1998 


August 9 


1 


47 


± 





01 


BD46AG 


1998 


August 23 


1 


49 


± 





01 


BD46AH 


1998 


September 9 


1 


52 


± 





01 


BD46AI 


1998 


September 25 


1 


55 


± 





01 


BD46AJ 


1998 


October 14 


1 


59 


± 





01 


BD46AK 


1998 


October 29 


1 


61 


± 





01 


BD46AL 


1998 


November 17 


1 


65 


± 





01 


BD46AM 


1998 


December 6 


1 


68 


± 





01 


BD46AN 


1998 


December 23 


1 


71 


± 





01 


BD46AO 


1999 


January 5 


1 


74 


± 





01 


BD46AP 


1999 


January 23 


1 


77 


± 





01 


BD46AQ 


1999 


February 6 


1 


79 


± 





01 


BD46AR 


1999 


February 19 


1 


82 


± 





01 



'^The optical phase is computed using the 
optical maximum at MJD = 50773 cited by 



Gray. Humphreys, fc Yated (119991 ). and assuming 



their quoted uncertainty A(/) ~ 0.01. A rnean 



period 



of 557.4 days is adopted (IKholopov et al.lll985l ) 



^Failed epoch, that could not be successfully re- 
duced in full. 



'^This epoch was not scheduled. 
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Table 2. VLA observations of J0359+509 





Source 


Date 


Start 


End 


X 




Config.'^ 


Epoch 






(lAT) 


(lAT) 


(deg)^ 


{deg)b 




A 


J0359+509 


1997 May 24 


17:08:39 


20:01:20 


124.8 


5.5 


B 


C 


J0359+509 


1997 Jun 23 


17:37:39 


19:00:20 


112.1 


3.6 


CB 


E 


J0359+509 


1997 Jul 24 


12:35:00 


14:01:49 


130.2 


2.6 


C 


G 


J0359+509 


1997 Aug 16 


13:02:20 


14:31:09 


143.1 


2.4 


C 


J 


J0359+509 


1997 Sep 30 


10:38:20 


11:34:10 


157.1 


0.6 


CD 


K" 


J0359+509 


1997 Oct 13 


07:48:49 


09:03:19 


161.8 


3.1 


CD 


L 


J0359+509 


1997 Oct 31 


07:04:50 


08:01:59 


144.0 


0.8 


D 


T 


J0359+509 


1998 Feb 08 


00:02:00 


00:59:29 


111.6 


1.7 


AD 


wf 


J0359+509 


1998 Apr 06 


22:42:40 


00:40:59 






A 


AA 


J0359+509 


1998 May 28 


16:54:00 


18:17:00 


128.5 


8.4 


A 


AB 


J0359+509 


1998 Jun 08 


17:36:19 


19:33:30 


162.0 


24.0 


AB 


AE 


J0359+509 


1998 Jul 17 


15:31:39 


17:30:00 


144.4 


3.2 


B 


AF 


J0359+509 


1998 Aug 15 


15:37:10 


17:05:50 


149.3 


2.5 


B 


AI 


J0359+509 


1998 Sep 24 


11:01:10 


12:52:49 


150.0 


1.2 


B 


AJ 


J0359+509 


1998 Oct 20 


10:15:30 


11:46:20 


7.8 


18.2 


B 


AL 


J0359+509 


1998 Nov 15 


09:06:10 


11:03:49 


24.5 


8.6 


BC 


AN 


J0359+509 


1998 Dec 22 


07:38:59 


09:38:19 


75.6 


1.5 


C 


AP 


J0359+509 


1999 Jan 22 


00:11:00 


03:07:09 


60.1 


3.7 


C 


AQ 


J0359+509 


1999 Feb 07 


03:03:19 


04:33:50 


73.8 


4.7 


CD 


AR 


J0359+509 


1999 Mar 09 


01:07:39 


03:05:50 


75.2 


6.5 


D 



^Absolute electric vector position angle (EVPA) of the linearly-polarized emission from 
J0359+509, relative t o an assumed EVPA for 3C138 at 43 GHz of X3C138 = -14 deg 
jPerlev fc TavloiibOOSh . 

''Estimated standard error in Xi derived from independent anlyses of each 50 MHz VLA 
continuum spectral window, as -^Ixi ~ X2\- 

'^VLA configuration code. 

<^ Nearest BD46 epoch in Table [T] 

'^EVPA measurement from VLBA project BK50 (Mitra et al. , in preparation). 
'Data were not reducible at this VLA epoch. 
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Optical Light Curve of TX Cam 
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5.05x10^ 5.1x10' 
Modified Julian Day 
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Fig. 1. — The optical light curve of TX Cam as provided by the AAVSO. The markers at 
the base of the figure indicate the dates of the VLBA 43 GHz observations. Reproduced 
from Figure 1 of Paper I. 
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Fig. 2.— The measured VIA absolute EVPA of J0359+509 {upper panel), the VLBA EVPA 
{middle panel), and the difference {lower panel), equivalent to the EVPA rotation needed to 
calibrate the absolute EVPA of the VLBA polarization images. 
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Fi§. 3. For each epoch, a contour plot of the enhanced contrast Stokes / image at levels [1,2,5,10,20,40,80,100] % of the peak. Vectors 
are overlaid proportional to linearly-polarized intensity (on a scale 10 mas = 18.08 Jy/beam) and drawn at a position angle of the EVPA. All 
Stokes parameters (/, Q, U) are summed over velocity. The figures are 1360 X 1360 pixels in size with a pixel spacing of 50pas. An fixed elliptical 
Gaussian restoring beam of 540 X 420/Ltas at a position angle of 20 deg was used for all images, representative of the point-spread function for the 
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Fig. 5.— 



Same as Figure [31 
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Fig. 6.— 



Same as Figure [31 
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Fig. 8. — Same as Figure [31 
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Fig. 9.— 



Same as Figure [31 




Fig. 10. — The mean fractional linear polarization magnitude as a function of stellar pulsat 
phase. 



Fig. 11. — A time-series of expanded sub-images of the isolated SiO maser component visible 
in the far north of images between epochs {S-AH}, faUing in toward the northern shell. The 
plot parameters are the same as for Figure [3l except for a more finely-sampled set of enhanced 
contrast Stokes / contour levels at (^io~^+°-i'^ for k=l,30.) 
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Fig. 12. — A sequence of weighted-mean component positions for the northern component 
trajectory shown in Figure [TTl Error bars in x and y are shown for each weighted- mean 
position, and the EVPA is drawn at each epoch as a bold vector, of uniform, constant length 
over epoch. 
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Fig. 13. — A sequence of weighted-mean component positions for the northern-eastern com- 
ponent discussed in the text. Error bars in x and y are shown for each weighted-mean 
position, and the EVPA is drawn at each epoch as a bold vector, of uniform, constant length 
over epoch. 
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Fig. 14. — The best-fit linear projected trajectory to the northern component positions in 
Figure [121 extrapolated across the full diameter of the SiO maser shell at epoch {Z}, here 
chosen as the mid-point epoch of the infalling component trajectory. 
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Fig. 15. — A time-series of expanded sub-images of the isolated SiO maser component visible 
in the eastern part of images between epochs {AA-AL}, falling in toward the eastern shell. 
The plot parameters are the same as Figure [HI 
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Fig. 16. — A sequence of weighted-mean component positions for the eastern component 
trajectory shown in Figure [151 Error bars in x and y are shown for each weighted- mean 
position, and the EVPA is drawn at each epoch as a bold vector, of uniform, constant length 
over epoch. 



